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Abstract 

The need may arise to be able to simulate the migration of groundwater nanoparticles through the ground. 
Transportation velocities of nanoparticles are different from that of water and depend on many processes that occur 
during migration. Unstable nanoparticles, such as zero-valent iron nanoparticles, are especially slowed down by 
aggregation between them. The aggregation occurs when attracting forces outweigh repulsive forces between the 
particles. In the case of iron nanoparticles that are used for remediation, magnetic forces between particles contribute 
to attractive forces and nanoparticles aggregate rapidly. This paper describes the addition of attractive magnetic 
forces and repulsive electrostatic forces between particles (by 'particle', we mean both single nanoparticles and 
created aggregates) into a basic model of aggregation which is commonly used. This model is created on the basis of 
the flow of particles in the proximity of observed particles that gives the rate of aggregation of the observed particle. 
By using a limit distance that has been described in our previous work, the flow of particles around one particle is 
observed in larger spacing between the particles. Attractive magnetic forces between particles draw the particles into 
closer proximity and result in aggregation. This model fits more closely with rapid aggregation which occurs between 
magnetic nanoparticles. 
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Background 

There are a lot of types of nanoparticles and colloidal 
particles in groundwater [1], Some of them are formed 
naturally, others are generated synthetically and put into 
the ground by humans. Not only is the reactivity of 
particles important, but also their migration properties 
are examined. For example, natural bentonite colloids 
are released as a consequence of bentonite disposal of 
radioactive wastes and could carry adsorbed radionuclides 
in groundwater through granite [2,3]. Zero-valent iron 
nanoparticles are produced [4-6] and injected into the 
ground. Iron nanoparticles are able to migrate in ground- 
water through contaminated areas and remediate the pol- 
luted soils and water [7]. In the first case, the migration 
possibility is unwelcome. In the second case, the bet- 
ter the migration, the more effective of the remediation. 
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That is why a simulation of the migration of nanopar- 
ticles might be desirable. To simulate the migration of 
nanoparticles, the coefficient of transport retardation of 
the nanoparticles is needed. The coefficient represents the 
possible reduction in the rate of nanoparticle migration 
compared with nanoparticles with similar properties. The 
number of nanoparticles with similar properties changes 
over time due to aggregation and it influences the results 
of the migration experiments. A dynamic model of aggre- 
gation has to be included in the simulation programme 
of nanoparticle transport in flowing water. That is why 
mass transport coefficients are needed. The coefficients 
represent the frequency of nanoparticle collisions [8,9]. 

A commonly used model for mass transport coefficients 
[10,11] in describing aggregation is based on the colli- 
sions among nanoparticles caused by heat fluctuation, the 
velocity gradient of the water in which the nanoparticles 
are suspended and the different velocities of sedimenta- 
tion of nanoparticles of varying size. This model does 
not include the decrease in the rate of aggregation due 
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to repulsive electrostatic forces which occurs due to the 
electric double layer which builds up on nanoparticle sur- 
faces [12]. Further, in the case of magnetic nanoparticles, 
the aggregation rate is rapidly increased due to the attrac- 
tive magnetic forces between nanoparticles [4,13-16]. 
That is the reason why the model of aggregation has been 
expanded, enabling a more accurate model of aggrega- 
tion of iron nanoparticles in water to be achieved. The 
paper describes the extension of the mass transport coef- 
ficients by the attractive magnetic forces and repulsive 
electrostatic forces between the nanoparticles. 

Methods 

A model of nanoparticle aggregation 

Particles aggregate easily in groundwater. They create 
clumps of particles up to the size of several microme- 
tres [15] that cohere and reduce the ability of particles to 
migrate through the pores on the ground. The aggrega- 
tion of the particles is caused by processes that generally 
occur during particle migration. The reduction in mobil- 
ity can be formulated by a rate of aggregation given by 
mass transport coefficients P (m^s~^) [9,10]. The coeffi- 
cients give a probability Pij for the creation of an aggre- 
gate from particle / and particle ; with concentrations 
Hi, Hj of particles /, /, respectively (Equation 1). Parti- 
cle / means the aggregate is created from / elementary 
nanoparticles. 

i^ij = (1) 
Aj = A)+A-+Aj- (2) 

The coefficient (Equation 2) is given by the sum of mass 
transport coefficients of Brownian diffusion ^j, velocity 

gradient A? and sedimentation fify The concept is adopted 
from [10]. 

In the case of small nanoparticles, temperature fluctua- 
tion of particles has a significant effect on particle aggrega- 
tion [17]. Brownian diffusion causes a random movement 
of the particles and it facilitates aggregation. The mass 
transport coefficient for the Brownian diffusion [10] is 

1 2kBT (d^ + d^ 
^ 3ri didj 

where stands for Boltzmann constant, T denotes the 
absolute temperature, rj is the viscosity of the medium, 
and di is the diameter of the particle /. 

Another process causing aggregation is the drifting of 
nanoparticles in water. Water flowing through a pore of 
soil has a velocity profile. In the middle of the pore, the 
velocity of water is highest. Since the particles have differ- 
ent velocities, according to their location in the flow, the 
particles can move close together and create an aggregate. 



The mass transport coefficient for the velocity gradients 
of particles [10] is 

fil = ^G(di + djf, (4) 

where G is the average velocity gradient in a pore. 

Particles settle due to gravitational forces. The veloc- 
ity of the sedimentation varies for different aggregates 
depending on their size, so particles can move closer 
together and aggregate. The mass transport coefficient for 
the sedimentation [10] is 

A) = ^(Qp- Q) (di + d^f \df - df\, (5) 

where g is the acceleration due to gravity, q is the density 
of the medium, and ^p is the density of the aggregating 
particles. 

The magnetic properties of nanoparticles 

Because of the composition of nanoparticles, every 
nanoparticle has a non-zero vector of magnetization. 
According to [15], TOD A iron nanoparticles produced 
by the Japanese company Toda Kogyo Corp. (Hiroshima, 
Japan) [5], with diameter of 40 nm have saturation mag- 
netization 570 kA/m. This is the value for a substance 
composed of nanoparticles containing 14.3% of Fe^ and 
85.7% of Fe304. We use these data for our model. There- 
fore, we assume the same size magnetization vector for all 
nanoparticles. 

Our model of a magnetic field around an iron nanopar- 
ticle is based on the model of the magnetic field around a 
magnet described in [18]. The electromagnetic potential 
in the point r near a permanent magnet of volume V is 
equal to 

0(r) = j ^dV. (6) 

V 

where M is the magnetization vector at the point dV, the 
vector R is the difference between source of the magnetic 
field dV and the point r, R is the length of R. 

The intensity of the magnetic field H can be subse- 
quently computed as 

H(r) = -grad(0(r)). (7) 

Finally, the magnetic force between the source of the 
intensity of magnetic field H and a permanent magnet of 
volume V with a magnetization vector Mq at the point r is 
equal to 

F(r) = -j (Mo- grad)H(r)dy. (8) 

V 

In our previous work [19], the scalar potential of the 
magnetic field around one homogeneous spherical iron 



Rosicka and Sembera Nanoscale Research Letters 201 3, 8:20 
http://www.nanoscalereslett.eom/content/8/1/20 



Page 3 of 9 



nanoparticle with radius a located at the point (0, 0, 0) was 
derived as follows: 



0(r) =M 



27t 7t a 

III 

0 0 0 



(;V3 - /cos(^))/2 sin(6>) 



(9) 



where a is the radius of the nanoparticle, and (»%i,»%2>«^3) 
are the coordinates of the point r. Here, the direction of 
the magnetization vector M is set towards X'^^ and M is the 
magnitude of the vector M. 

From Equations 7 and 8, the analytical computation of 
the magnetic force between two iron nanoparticles can 
be obtained. Since nanoparticles aggregate, the magnetic 
force between aggregates must be derived. One aggregate 
can be composed of millions of nanoparticles. It would be 
time-consuming and very difficult to analytically compute 
all these forces. As a consequence, the forces are com- 
puted numerically, either as a sum of the magnetic forces 
between every nanoparticle in one aggregate with every 
nanoparticle in the second aggregate 



«2 



F = ^ V{M2] • grad)grad0(r2j), 



(10) 



or as one magnetic force between two averaged aggregates 
[20]. 



F = K2(M2A • grad)grad(/)(R,MiA, ^^a). 



(11) 



where V = ^na^ is the volume of a nanoparticle, r2j is 
the location of the centre of the ;-th nanoparticle in the 
second aggregate, M2j is the magnetization vector of the 
/-th nanoparticle in the second aggregate, Mia and M2A 
are the averaged magnetization vectors (Equation 12) of 
the first and the second aggregate respectively, and V2 = 
^/n2a is the volume of the second aggregate. 

The averaged aggregate is a big homogeneous particle 
with its direction of magnetization vectors Ma which is 
computed as a vector sum of the magnetization vectors 
of all nanoparticles in the aggregate Ma and computed as 
an average of the sizes of all nanoparticles divided by the 
number of nanoparticles in the aggregate n. 



Ma 



(12) 



The structure of aggregates 

When particles aggregate due to magnetic forces, the rate 
of aggregation depends on the magnetization vectors of 
the aggregating particles and on the distance between the 
particles. The rate of aggregation changes with the chang- 
ing number of nanoparticles within the aggregates, that is, 
the changing scale of the structure by order. The model 
which has been chosen for the structure of an aggregate 
is a sphere with randomly located nanoparticles within 



the aggregate, either with random directions of magne- 
tization vectors for every nanoparticle; or with the same 
direction of magnetization vectors for all nanoparticles in 
the aggregate. Aggregate structures were assessed in pre- 
vious work [21]. A more accurate assessment of the most 
probable structure of an aggregate was performed for this 
paper in section 'The structure of an aggregate based on 
interaction energy! 

The electrostatic properties of nanoparticles 

In an electrolyte, a surface charge builds up on the 
nanoparticle surface. The surface charge depends on its 
zeta potential (see e.g. [22]) which is measurable. The 
zeta potential strongly depends on the pH of the water. 
The results of this dependence were measured using the 
Malvern ZetaSizer (Malvern Instruments Inc, Malvern, 
Worcestershire, UK) as published in [19]. From the zeta 
potential, the surface potential can be computed, based on 
the electrical double layer [23,24] 



a = —^8€o€r:cRg Tsiuh 



2Rf 



(13) 



where a is the surface charge density of the particle, c is 
the molar electrolyte concentration, Rg is the molar gas 
constant, F is Faraday s constant, Z is the charge number 
and f is the electrostatic potential. The electrostatic force 
between two particles is equal to 



46o6r D2 



(14) 



where D is the distance between the particles / and 
The electrostatic forces repel nanoparticles with the same 
polarity and cause a reduction in the rate of aggregation. 
Inclusion of the dependence is done in section 'The inclu- 
sion of the limit distance into mass transport coefficients! 

The limit distance 

The effect of magnetic forces on the rate of aggregation 
was assessed by one parameter - the limit distance Ld- 
This dimension expresses the range of magnetic forces 
between particles. The definition of this parameter is as 
follows: this is the distance from centre of an aggregate up 
to which attractive magnetic forces cause the aggregation 
between the aggregate and a particle placed in this range. 
Hence, in a range larger than the limit distance, other 
forces outweigh the magnetic forces (Figure 1). The limit 
distance Lj) can be defined as the distance of the point in 
which gravitation Fg and magnetic forces Fmg effecting on 
the aggregate are equal 



(15) 
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Figure 1 Sketch of the limit distance. A comparison of the forces 
acting on aggregates depicted by a two-dimensional figure. Inside 
the circle with diameter equal to the limit distance, the magnetic 
forces outweigh the gravitational force and aggregation occurs. 
Outside this, the aggregates settle. 



The limit distance takes the form 



(16) 



The magnetic force between two single domain mag- 
netic nanoparticles falls by the power of 4. In the case of 
aggregates, the fall depends on the structure of the aggre- 
gates and iteration of limit distance computation is needed 
[20]. 



^D,l = J ^D,0. 



(17) 



When including electrostatic forces, we define the limit 
distance as the distance where the repulsive magnetic 
forces is equal to the sum of attractive forces Fmg and Fc« 
As the effect of electrostatic forces falls by the power of 2, 



electrostatic forces can be included into the equilibrium 
of force in the following way [20] 



Fl(Ro) + 4F^F^^(Ro)-Fc(Ro) 
2K 



Ro> (18) 



Fl(Lj),o) + 4FgFmg(^D,o) - Fc(Lj),o) 



2F^ 



(19) 



The values of magnetization vector and surface charge 
were selected as follows: M = 570 kA/m; a = 2.5 x 
10~^ C/m^. We used these selected values for all the com- 
putations of the interaction energies and mass transport 
coefficients. 

Simulation software 

All the computations of magnetic forces, limit distance, 
electrostatic forces and mass transport coefficients were 
performed using Matlab R2009a software (MathWorks 
Inc, Natick, MA, USA). The computation was carried out 
for different sizes of aggregates / and mostly varying in 
the order of the number of nanoparticles that the aggre- 
gates were composed of. The magnetic forces between 
two aggregates were computed either by summation of 
the magnetic force between every nanoparticle in the first 
aggregate and every nanoparticle in the second aggregate 
(when the ratio Lj^/Rq expresses distance between the 
aggregates was lower than 15 [20]), or by the averaging 
of the first and second aggregates. Values for the mag- 
netization vector and surface charge were selected in the 
following way: M = 570 kA/m; a = 2.5 x 10"^ C/m^. For 
the velocity gradient, we chose the dimensionless value 50. 
We used these selected values for all the computations of 
the interaction energies and mass transport coefficients. 

Results and discussion 

The structure of an aggregate based on interaction energy 

To assess the most probable structures of aggregates, 
one can compute an interaction energy E between the 
nanoparticles which make up the aggregate, according to 
[25] 



-m B. 



(20) 



This is the potential energy of the magnetic moment 
m in the externally produced magnetic field B. Again, we 
assume the same magnetization vectors for all nanopar- 
ticles in the aggregates with value 570 kA/m [15]. Posi- 
tive interaction energy means repulsion of the magnetic 
moment from the magnetic field of another magnetic 
moment; negative interaction energy means attraction of 
the dipoles. By summation of the interaction energies 
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Figure 2 The interaction energies of different structures of aggregates. A comparison of the interaction energies of different structures of 
aggregates expressing tine rate of probability of tine structures (tine larger tine negative energy, the bigger the probability of structure). 



between every two nanoparticles in an aggregate, one can 
deduct the probability of stability of the different struc- 
tures of the aggregates (the higher the negative interaction 
energy, the higher the probability of the structure of the 
aggregate). 

The results of interaction energies are shown in Figure 2. 
The computed interaction energies are displayed for dif- 
ferent structures of aggregates (according to the schemes: 
Figures 3, 4, 5, 6). The Figure 2 is shown using a log- 
arithmic scale. The exact values of interaction energies 
for different structures of aggregate (Figures 3, 4, 5, 6) 
and the different numbers of nanoparticles making up 
the aggregates are in Table 1. Not the absolute values 
but the comparison between the values of the different 
structures is relevant. According to Figure 2, the most 
probable structure of aggregates for the small aggregates 
are chains and for the bigger aggregates, spherical clusters 
with the same direction of magnetization vectors of the 
nanoparticles which make up the aggregate. 

In their research, Phenrat et al. [15], aggregates of 
nanoscale zero-valent iron particles were measured using 




Figure 3 Diagram of a chain structure. A diagram of the chain 
structure of nanoparticles within an aggregate with schematic 
directions of the magnetization vectors of the nanoparticles. 



dynamic light scattering, optical microscopy and sed- 
imentation measurements. According to their results, 
firstly, the nanoparticles created clusters and subse- 
quently, these aggregates assemble themselves into fractal, 
chain-like clusters. We presume that it was because of the 
high concentration of nanoparticles that they used, and 
the very fast aggregation, first into chains and then into 
clusters, which lead to the measurement of only larger 
clusters in [15]. Our presumption that with larger num- 
bers of nanoparticles, spherical cluster is created which 




Figure 4 Diagram of a circular structure. A diagram of a circular 
structure of nanoparticles within an aggregate with schematic 
directions of the magnetization vectors of the nanoparticles. 
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Figure 5 Diagram of spherical structure. A diagram of a spherical structure of nanoparticles witliin an aggregate witli scliematic directions of tine 
magnetization vectors of tlie nanoparticles. 



leads to the supposition that at very high concentrations 
of particles, spherically structured aggregates only attach 
to each other, without changing their structure. This cor- 
responds to the observations of Phenrat et al. [15]: in high 
concentrations, first nanoparticles aggregate into clusters, 
then the created clusters aggregate into pairs or triplets, 
and finally into chain-like fractal aggregates. 

The inclusion of the limit distance into mass transport 
coefficients 

The basic model of aggregation as given in the section, 
'A model of nanoparticle aggregation', indicates the rate of 
aggregation caused by the collision of particles (in prox- 
imity, attractive forces outweigh the repulsive ones). We 
established a limit distance in which attractive forces out- 
weigh the repulsive ones. The magnetic forces attract 
particles closer to each other and then they aggregate due 
to attractive van der Waals forces. 

Mass transport coefficients (in Equations 3, 4, and 5) 
were derived on the basis of the flux of nanoparticles 
through an observed volume or circular area around a par- 
ticle. The area had a radius equal to sum of the radii of 
both particles. That means that the particles collide and 
aggregate. According to our supposition, the particles do 
not have to be in proximity to aggregate when attractive 



magnetic forces are acting between them. Therefore, the 
mass transport coefficients are computed as flux through 
the spherical or circular area around a particle with a 
diameter equal to the limit distance: 

i,mg _ 4/:b T ( 1 
Pij - 



3^ 



+ 



d) 



(21) 



Q2,mg 
Q3,mg 



GL 



18 r] 



(22) 



(23) 



where pf^'^^, pfj^^^ and pf^'^^, stand for the mass trans- 
port coefficient of Brownian motion, the velocity gradient, 
and sedimentation respectively, with the inclusion of mag- 
netic forces between particles. The results of this change 
in mass transport coefficients are discussed in the next 
section - 'A comparison of the rate of aggregation with and 
without the effect of electrostatic and magnetic forces! 

A comparison of the rate of aggregation with and without the 
effect of electrostatic and magnetic forces 

The comparison was carried out using an extreme case 
with a spherical aggregate structure with the same direc- 
tion of magnetization vectors of all nanoparticles within 
the aggregates. The aggregation is highest in this case 
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because attractive magnetic forces attract the aggre- 
gates and the rate of aggregation is significantly higher 
(Figure 7). Table 2 contains a comparison of mass trans- 
port coefficients computed by primary model, mass 
transport coefficients computed in distance Lj) including 



magnetic forces and mass transport coefficients computed 
in distance Lj) including both magnetic and electrostatic 
forces. The computation of Id was performed by aver- 
aging the magnetic forces for particles with ratio Ld/^o 
higher than 15; otherwise, the computation of magnetic 



Table 1 Interaction energies of different structures of 
aggregates 



Number of nanoparticles [1 ] 


Structure 


Energy/ fi (eV) 


2 


Chain 


273 


3 


Chain 


588 


8 


Cube 


903 


8 


Sphere 


1,449 


8 


Circle 


2,184 


8 


Chain 


2,688 


27 


Chain 


3, 780 


27 


Sphere 


8,400 


29 


Cube 


8,400 


343 


Cube 


56,700 


343 


Chain 


1 09, 200 


343 


Sphere 


1 84, 800 



Computed interaction energies divided by the permittivity constant for different 
structures of aggregates (according to the diagrams in Figures 3, 4 ,5 ,6) and for 
different numbers of nanoparticles within the aggregates. 



Number of nanoparticles in aggregate [1] 

100 1000 10000 100000 1000000 



lE-14 
J [mVs] 



V 

MTC 



-MTC with magnetic forces 



MTC with magnetic and 
electrostatic forces 



Figure 7 Mass transport coefficients (MTC) comparison. A 

comparison of mass transport coefficients computed by the primary 
model, mass transport coefficients computed in distance 
including magnetic forces, and mass transport coefficients computed 
in distance including both magnetic forces and electrostatic 
forces. The MTC represents the sum of MTCs for Brownian motion, 
velocity gradient, and sedimentation. 
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Table 2 Comparison of mass transport coefficients 



/[1] 










^Dmg r -1 

/?o ^ 'J 


1 


1 


1.1 X 10-^^ 


3.1 X 10-^5 


2.9 X 10-^5 


78.9 


1 


10 


1.3 X 10-^7 


2.9 X 10-^5 


2.8 X 10-^5 


50.6 


1 


100 


1.9 X 10-^7 


2.8 X 10-^5 


2.7 X 10-^5 


28.4 


1 


1,000 


3.4 X 10-^7 


2.7 X 10-^^ 


2.7 X 10-^^ 


14.6 


1 


10,000 


7.3 X 10-^7 


2.8 X 10-^5 


2.8 X 10-^5 


7.1 


1 


100,000 


2.2 X 10-^6 


3.1 X 10-^5 


3.0 X 10-^5 


3.4 


1 


1,000,000 


1.4 X 10-^5 


4.2 X 10-^5 


4.2 X 10-^5 


1.6 


10 


10 


1.1 X 10-^7 


1.4 X 10-^4 


1.3 X 10-^4 


65.6 


10 


100 


1.3 X 10-^7 


1.3 X 10-^4 


1.3 X 10-^4 


42.0 


10 


1,000 


2.0 X 10-^7 


1.3 X 10-^4 


1.3 X 10-^4 


23.5 


10 


10,000 


4.2 X 10-^^ 


1.3 X 10-^4 


1.3 X 10-^4 


12.1 


10 


100,000 


1.6 X 10-^^ 


6.9 X 10-^4 


6.8 X 10-^4 


10.2 


10 


1,000,000 


1.3 X 10-^5 


2.5 X 10-^4 


2.5 X 10-^4 


3.2 


100 


100 


1.2 X 10-^7 


7.1 X 10-^4 


6.9 X 10-^4 


54.4 


100 


1,000 


1.5 X 10-^7 


7.1 X 10-^4 


7.0 X 10-^4 


34.7 


100 


1 0, 000 


3.0 X 10-^7 


7.2 X 10-^4 


7.1 X 10-^4 


19.4 


100 


100,000 


1.4 X 10-^^ 


7.0 X 10-^3 


7.0 X 10-^3 


21.1 


100 


1,000,000 


1.3 X 10-^5 


1.9 X 10-^3 


1.9 X 10-^3 


6.4 


1,000 


1,000 


1.5 X 10-^7 


4.0 X 10-^3 


3.9 X 10-^3 


45.1 


1,000 


10,000 


3.2 X 10-^7 


4.0 X 10-^3 


4.0 X 10-^3 


28.7 


1,000 


100,000 


1.5 X 10-^^ 


4.1 X 10-^3 


4.1 X 10-^3 


16.1 


1,000 


1,000,000 


1.4 X 10-^5 


1.3 X 10-^2 


1.3 X 10-^2 


11.8 


10,000 


10,000 


5.4 X 10-^7 


2.2 X 10-^2 


2.2 X 10-^2 


37.3 


10,000 


100,000 


2.2 X 10-^6 


2.3 X 10-^2 


2.3 X 10-^2 


23.7 


1 0, 000 


1,000, 000 


1.8 X 10-^5 


2.4 X 10-^2 


2.4 X 10-^2 


13.3 


100, 000 


100,000 


4.4 X 10-^^ 


1.3 X 10-" 


1.3 X 10-" 


30.8 


100,000 


1,000,000 


2.7 X 10-^5 


1.3 X 10-" 


1.3 X 10-" 


19.6 



A comparison of mass transport coefficients computed by the primary model fi, mass transport coefficients computed in distance Lq including magnetic forces fimg, 
and mass transport coefficients computed in distance including both magnetic forces and electrostatic forces ^6^5. The represents the sum of the mass transport 
coefficients for Brownian motion, velocity gradient and sedimentation. Computation of was performed by averaging of the magnetic forces (11) for particles with 
ratio /-D /Ro higher than 1 5; otherwise, the computation of magnetic forces was performed accurately by summation (1 0) (for more information, see [20]). Results were 
computed using the following values: M = 570 kA/m; a = 2.5 x 1 0"^ C/m^; and G = 50. 



forces was done accurately by summation (for more infor- 
mation see [20]). The values in Table 2 are computed with 
values M = 570 kA/m; a = 2.5 - 10'^ C/m^; G = 50. 
According to the results in Table 2 for the chosen values 
of variables, the attractive magnetic forces between iron 
nanoparticles have a large effect on the rate of aggrega- 
tion. The mass transport coefficients are much higher and 
the aggregation probability increases, which corresponds 
to our expectations. 

Discussion 

In future work, the system of grouping of particles accord- 
ing to their size will be derived for the new extended 
mass transport coefficients including electrostatic and 



magnetic forces. The groups will represent particles with 
similar transport properties (small particles are easily 
transportable, large particles remain in the pores in the 
ground) and a model of aggregation over time will be 
developed. The model will be compared with the mea- 
suring of aggregation of zero-valent iron nanoparticles 
in time. 

Subsequently, the limit distance should be derived for 
the equilibrium of all forces acting on particles depending 
on specific conditions. 

Conclusions 

In the case of magnetic nanoparticles with non-zero sur- 
face charges migrating through the ground, a basic model 
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of interaction between nanoparticles described by the 
probability of collision due to Brownian motion, velocity 
gradient, and sedimentation is insufficient. In our pre- 
vious work, we derived the level of effect of repulsive 
electrostatic forces between the nanoparticles, and we 
assessed the level of effect of the attractive magnetic forces 
between magnetic nanoparticles. In this paper, we sum- 
marised the findings and included it into an analytical 
model of collisions between magnetic nanoparticles. Due 
to attractive magnetic forces, the rate of aggregation is sig- 
nificantly higher, whereas the repulsive electrostatic forces 
are almost negligible. One can suppose that with other 
realistic selections of values of magnetization vector or 
surface charge, this trend would not change dramatically. 
This modified model of aggregation can better explain 
the rapid aggregation of zero-valent iron nanoparticles 
that is observed. This can help with the simulation of the 
migration of undissolved particles in groundwater. 
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